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Abstract 

The SuperCDMS experiment implements transition-edge sensors to measure athermal phonons produced 
by nuclear recoils in germanium crystals. We discuss a numerical simulation of these TES devices and a 
procedure for tuning the free model parameters to data, which reproduces superconducting-to-normal phase 
separation within a TES. This tuning provides insight into the behavior of the TESs, allows us to study the 
phase-separation length to optimize our detector design, and is integrated into a more complete simulation 
of the phonon and charge physics of SuperCDMS detectors. 



1. Introduction 

Transition-edge sensors (TESs) function as extremely sensitive calorimeters and are used in applications 
such as dark matter detection (TJ |2] , x-ray astrophysics [31 13] , and measurements of the cosmic microwave 
background [3 [S] . A TES consists of a superconducting material typically held under a constant voltage 
bias that maintains it in the superconducting transition. Negative electrothermal feedback stabilizes the 
sensor against small perturbations in current or temperature [7 . In this arrangement, a TES is a stable 
device that experiences a measurable current change in response to an energy input. Applications of TESs 
exploit this behavior to achieve high energy sensitivity to small signals [5]. 

One technical complication that can arise when operating TESs is the occurrance of superconducting-to- 
normal phase separation within the device [9]. Instead of being homogeneous, different parts of a TES may 
have different temperatures or thermodynamic phases. This typically occurs when there is weak thermal 
coupling between different regions of the TES. For example, the ends of a TES whose length is much greater 
than its width may be in different phases. 

To heuristically understand why phase separation occurs, consider a long, voltage-biased TES, as shown in 
Figure[lJ We can consider it as consisting of two thermally decoupled regions, denoted A and B. If a thermal 
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Figure 1 : If length of a TES is much longer than its width, then thermal diffusion between the two ends of the device can be very 
small. In this situation, the TES behaves similar to two thermally disconnected TESs in series, and superconducting-normal 
phase-separation results. 



fluctuation causes the temperature of A to rise, then tlie resistance will also increase and, because the TES is 
voltage biased, the current will decrease. This decreases the Joule heating and subsequently the temperature 
in B. Region A thus tends toward the normal phase and region B tends toward the superconducting phase. 
Thermal diffusion along the length of the TES counteracts this effect, but diffusion may be insufficient to 
prevent phase separation in a long TES. 

Phase separation has several negative effects in applications. The current response to a voltage or 
temperature impulse is different for phase-separated and phase-uniform cases. Phase-separated TESs also 
have higher noise and reduced sensitivity |10] . Internal thermal fluctuations become more significant because 
of the larger temperature gradient across the device. These internal thermal fluctuations produce additional 
noise in TESs [II] 

The SuperCDMS experiment searches for evidence of weakly-interacting massive particles (WIMPs) by 
detecting nuclear recoils in high-purity germanium detectors ^12, . Nuclear recoils produce phonons and break 
electron-hole pairs in the Ge. The phonons are detected at the top and bottom faces of the Ge crystals using 
tungsten-aluminum quasiparticle-trap-assisted electrothermal feedback transition-edge sensors (QETs) [13] • 
The phonons from the Ge substrate are collected by superconducting Al "fins" that cover a substantial 
fraction of the crystal surface. The absorbed phonons break Cooper pairs and produce quasiparticles. 
Quasiparticles diffuse through the Al to the interface with the W TES, where W and Al overlap to form a 
bilayer, as shown in Figure [2] The quasiparticle energy in the W is lower than the energy in the bilayer, 
which is in turn lower than in the Al. As a result, the quasiparticles are ultimately trapped and release 
their energy in the bilayer region and the W TES. The energy from the quasiparticles is then detected as a 
change in voltage across the TES due to an increase in temperature ^14) . A large number of TESs are read 
out in parallel in each of the eight channels of a detector. 

An important point is that some of the phonons in the Ge substrate, so-called athermal phonons, carry the 
position information of the initial event. Thus, maximizing the detected signal from this phonon population 
is an important design goal. Doing so requires maximizing the phonon collection area on the surface and 
maintaining sufficient bandwidth to resolve the signal. Increasing the phonon collection area is achieved 
simply by increasing the number of TESs and their attached Al fins. The bandwidth, on the other hand, is 
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Figure 2: Once quasiparticles are produced in tlie Al, tliey diffuse into the overlap bilayer region where the gap energy is 
smaller. The quasiparticles then decay by releasing phonons, and they become trapped in the overlap region. The same process 
occurs when the quasiparticles diffuse into W. 



set by the L/R time-constant that determines the frequency at which the detector response roUs off. Since 
the inductance L is fixed by the SQUIDs that are used to read out the TESs, the resistance of the TESs is 
set in the design to match the desired bandwidth. In SuperCDMS detectors, the actual phonon channels 
each consist of many hundred TESs arranged in parallel. In order to maintain the desired channel resistance 
as more TESs are added, the length of each TES must be increased. Sufficiently long TESs can separate 
into two phases, degrading the signal to noise. Thus, optimizing the detectors involves a trade-off between 
maximizing phonon collection and minimizing the phase separation that can result from using long TESs. 

In this paper, we discuss a simple simulation of the SuperCDMS TESs which allows us to study phase 
separation in the detectors. The model has several parameters: some that can be measured directly and 
others that must be tuned using diagnostic data. We describe the effects of each model parameter on the 
simulation, and describe a procedure for tuning the simulation to better understand and predict detector 
performance. 

2. Simulation 

To simulate a TES that could be phase-separated, we model a single TES as a one-dimensional array of 
nodes connected in series [9|. The length of a TES in a SuperCDMS detector is the distance between the 
first and last nodes of this one-dimensional array. We let each node have its own temperature and resistance, 
so that the two ends of the TES can have different phases. 

In the model and experiments, each TES can be characterized by several static and tunable parameters, 
as discussed below |15) . 

1. Parasitic resistance (Rp) - In the superconducting phase, the circuit has a residual parasitic resistance 
due to impurities and resistance in the leads that read out the TESs. 
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Figure 3: Diagram of the TES circuit (a) and the Thevenin equivalent diagram (b). From [8]. 

2. Normal resistivity {pn) - Each W TES has a characteristic normal resistance when it is just above the 
critical temperature, although small fabrication differences may cause this to be different for different 
TESs. We denote the normal resistance of a TES by 

3. Electron-phonon coupling parameter (S) - This parameterizes the thermal coupling between the elec- 
tron system and the phonon system within the W TES and W spurs that connect to the Al fins. It has 
been measured to be approximately 3.2 x 10^ Wm~'^K~^, in the absence of the Al fins that connect to 
the Ge substrate and the TES [16]. The geometry of these connections can vary somewhat between 
TESs, and the effect of this variation and the volume of Al material on E of the TES is difficult to 
quantify a priori. We therefore allow this parameter to be tuned in order to match the simulation to 
diagnostic data. 

4. Specific heat capacity (C) - The specific heat capacity of a TES used in SuperCDMS detectors ideally 
should be the value for W. However, the presence of the Al fins changes the effective specific heat of 
the device in ways that are difficult to quantify a priori. As a result, we also allow this parameter to 
vary. 

5. Diffusion constant (D) - This parameterizes the one-dimensional diffusion along the length of a TES. 
Because of possible sub-micron-scale variations in the photolithographically patterned W TESs on the 
detectors, we treat this as a tunable parameter. 

6. Critical temperature (Tc) - The critical temperature of a TES is the temperature of the superconducting 
transition, which may vary slightly from one TES to another. This can be measured to within about 
1 mK. The critical current of a TES (the current needed to drive the TES normal) is measured at 
several different temperatures and the critical current vs. temperature data is extrapolated to zero 
current. The temperature corresponding to zero current is taken to be the critical temperature of the 
TES. 

7. Transition width (ATc) - This describes the width of the superconducting transition, and it is related 
to the 10-90% transition width by ATio-go = 2Arc tanh"\0.8). For the simulation, we assume 
that ATio_go — 0.35 niK, consistent with measurements. This value is determined from resistance 
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measurements. 

8. Spatial distribution of critical temperatures (characterized by variance a'^) - In the model, the Tc of 
each TES is ahowed to vary across the surface of the detector. Due to the process by which the W is 
patterned on the Ge, small variations in the value of Tc across the detector can occur. In particular, 
small but nonzero radial gradients and linear gradients across the detector surfaces have been measured 
in some detectors. Measurements of the detectors considered here indicated a small random variation 
of Tc between individual TESs, but no systematic gradient. For the purposes of this discussion, we 
therefore define the Tc of each TES using a Gaussian with mean at the measured Tc and a tunable 
variance a^. 

In our model, heat can be transferred to and from an individual TES by two processes: heat transfer to 
the substrate and Joule heating. Radiative effects are negligible. In order to model the heat transfer to the 
substrate, the functional form of the resistance of the TES as a function of temperature and current must 
be assumed. We take Ri{Ti,Ii) to have the following form, motivated by Ginzburg-Landau theory near Tc 



Ri{Ti,Ii) — 



l + tanh(^^M_^ + (j(i)/j„)2/3^ 



(1) 



where /q is a constant, and the subscript i labels the nodes into which each TES is subdivided. Ti(i) is the 
temperature of the electron system of the ith node at time and Ii{t) is the current through the ith node 
at time t. Joule heating and cooling of the TES to the substrate then gives the net power transfer to the 
TES, 

= VY.{Tl,^-TSf)+I^R,{t), 

where C and V are the specific heat capacity and specific volume of W for a single TES node, Tf,ath is the 
temperature of the constant temperature of the phonon system. At is the time step of the simulation, and 
S is the electron-phonon coupling parameter described above. The first term gives the cooling to the Ge 
substrate and is motivated by [16] . The second term gives the Joule heating. 

Additionally, heat can be transferred between neighboring nodes within one TES by diffusion. This 
process is described by a discrete ID diffusion equation 

^ r,+i(0-2r,(t) + r,_ift) _ r,(t)-r,(t-i) 

^ (Ax)2 - At ' 

where D is the diffusion constant, and Ax is the width of each node. 

Each SuperCDMS detector has eight independent phonon channels (four on the top surface and four 
on the bottom surface), each of which contains 500 TESs in parallel and wired in a constant-voltage 
configuration. Real data obtained with a SuperCDMS detector therefore refiects the ensemble response of 
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Figure 4: Current response in the TES due to a slowly decreasing bias voltage in the phase-uniform case. 

many different TESs. However, given the uniformity of the patterned TESs and the quahty of the data 
observed, the paraUel array has been shown to behave quahtatively the same as a single TES. 

3. Tuning with Ib-Is curves 

To use the simulation that we developed, the parameters described in the previous section must be 
tuned to data. A bias current-sensor current {Ii,-Is) curve is one useful test that we perform in the lab, 
which allows us to tune Rp, i?„, S, and a^. This test involves applying a linearly decreasing bias current 
to the TES circuit and measuring the current response. Note that the TES itself is voltage biased, but 
we actually monitor the bias current through the parallel TES circuit (see Figure |3]). The It-Is curve is 
therefore equivalent to an I-V curve. The rate of decrease in current is chosen to be slow compared to all 
time scales of the detector response so that the TES is approximately in thermal equilibrium with the Ge 
substrate at all times during the test. Since the device is in thermal equilibrium, we can recast the lb-Is 
data into a plot of power dissipation as a function of applied bias current, shown in Figure [Sj 

At high applied bias current (e.g. 0.3 mA), a channel of TESs has an approximately constant resistance 
from the parasitic resistance (~15mri) and the normal resistance of the W (^0.7fi), and it obeys Ohm's 
law. Similarly, at low current (e.g. 0.02 mA), the TES obeys Ohm's law and has an approximately constant 
resistance given by the parasitic resistance alone. This is evident from the linear regions of Figure |4] As the 
bias current decreases from high current and the TES begins to enter the transition, the TES temperature 
is stabilized by electrothermal feedback and slowly progresses through the transition. As the temperature 
decreases, the resistance decreases, and the current through the TES increases, producing the spike in 
Figure |4] Since the TES is in equilibrium and its temperature is nearly constant in the transition, the power 
dissipated to the substrate is constant, producing the plateau seen in Figure [5| 
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Figure 5: Power dissipated through the TES as a function of a slowly decreasing bias voltage in the phase-uniform case. 
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Figure 6: Experimental results (solid blue) and simulation (dashed red) for the current response of a phase-separated TES as 
a function of the TES bias current. 



In the phase-separated case, part of the TES is the normal phase and another part is the superconducting 
phase. When this occurs, the shape of the Ib-Ig curve changes in the transition, as shown in Figure [6j The 
plateau region of the power curve becomes considerably more sloped, as shown in Figure [7] At some point 
in the transition the TES switches from a phase-uniform to phase-separated state. In the phase-separated 
state, a significant fraction of the TES is at fixed resistance-either the normal resistance or nearly zero 
resistance. As the bias current decreases, the power dissipation also decreases. 

To tune Rp, Rn, S, and ct^ with these data, we attempt to find the parameters that best match the 
bias current versus power curve for simulation and data (e.g. Figure 4). First, the parasitic and normal 
resistance are determined by performing a linear fit to the normal and superconducting regions of the Ib-Is 
curve and extracting the slope of each, as shown in Figures [4] and [6] for example. The parameters T, and 
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Figure 7: Experimental results (solid blue) and simulation (dashed red) for the power dissipated by a phase-separated TES as 
a function of TES bias current. 

are tuned by running the simulation for a matrix of parameter values. The best- fit values of the parameters 
are selected by the least-squares difference between the data and simulated It-Is curves in the transition 
region. This process is iterated to obtain the optimal parameters. 

Graphically, increasing S or Tc has the effect of raising the flat region of the plot of bias current versus 
power dissipated. This occurs because the heat transfer to the substrate is larger and the transition is 
reached at a higher bias current. 

The data from Ib-Is curves alone cannot break all of the degeneracies between the model parameters. 
There is considerable degeneracy between S and Tc, for example. Additionally, C and D have only a small 
effect on the shape of the lb-Is curve, and this method is consequently fairly insensitive to these parameters. 
For this reason, we take the approximate from critical current measurements, and we do not tune C or 
D with the lb-Is curves. 

While the model fit to lb-Is curves is generally good, it should be noted that systematic error is introduced 
through our assumption that the form of the resistance curve in the transition is given by equation ([T]) . This 
is a phenomenological parameterization that produces acceptable results, but it does not represent the true 
form of the resistance curve. Differences between this form and the true resistance curve contribute to the 
remaining discrepancy between simulated and measured curves. 

4. Tuning with sudden impulse 

To better tune the diffusion constant and the specific heat capacity, we study the detector response to 
a square wave current. We apply a square wave of 0.025 mA peak-to-peak amplitude on top of a constant 
bias current of 1.54 mA. A characteristic response associated with the rising edge of a square wave is 
shown in Figure m The sudden increase in the voltage across the TES initially increases the current with a 
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Figure 8: Experimental (solid blue) and simulated (dashed red) TES response to a square wave bias current. 

characteristic time scale L/R. At slower time scales, the TES temperature and its resistance increase through 
Joule heating, which decreases the current flow. Eventually a new quasi-equilibrium point is reached with 
both smaller current and larger resistance than originally. 

The specific heat controls the temperature change needed to restore internal equilibrium and the rate at 
which energy is exchanged to the substrate, while the diffusion constant affects the time to reach internal 
thermal equilibrium along the TES. For these reasons, the specific heat and diffusion constants strongly 
affect the decay time of the TES response to a square wave, and both parameters can be tuned by fitting 
data and simulation. The parameter S also has a significant effect on the decay time of the current response 
by controlling the rate at which the TES reaches thermal equilibrium with the substrate. 

We find, in particular, that increasing the diffusion constant tends to decrease the decay time. Increasing 
the specific heat, on the other hand, increases the decay time. A larger specific heat implies that the 
temperature excursions due to the changing bias current are smaller; therefore, the cooling power will be 
smaller and the time to return to equilibrium will be larger. The difference between the sensor current 
before the square wave edge and at a long time after is determined by the overall constant bias current 
applied: being biased at different points in the superconducting transition will result in different changes in 
the resistance of the TES when a square wave bias is applied. 

The response of a phase-uniform TES is characterized by a single exponential risetime and a single expo- 
nential falltime 8 . In contrast, the response of a phase-separated TESs is best described by a single expo- 
nential risetime and two exponential falltimes. The qualitative behavior of these pulses is well-reproduced 
by the simulation, as seen in Figure [8] although it is difficult to precisely match the parameters of the 
exponentials describing the data and simulated pulses. 
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Parameter 


characteristic value 


mean Tc 


90 mK 




0.35 mK 


(Tj' 


1 mK 


E 


4.8 X 10^ W K-5 


Pn 


1.2 X lO"'' n m 


D 


4.03 X 10-" m-2 s 


Cp 


37.0 J K-i m-3 



Table 1: Characteristic model parameters found from matching of /(,-/s curves and square-wave response. 



5. Phase-Separation Length 

The phase-separation length for a TES is the maximum length of a TES before it becomes phase- 
separated. There is a simple analytic expression for this length, assuming the Wiedemann- Franz law jl6j . 

where L = 2.44 x 10^^ W ft K^^ is the Lorentz constant, p„ is the normal resistivity, and a = OlogR/dlogT. 
This expression is not exact because Wiedemann-Franz holds only approximately as these temperatures. For 
the functional form of the resistance function that we assume in equation ([T]), we have 

T 1 

a(u) = —— ^ , (4) 

^ ' ATe(cosh2M)(l+tanhu) ^ ' 



where 



T~Tr, 



"=^^ + (//^)'/^ (5) 

and T and / are the temperature and current at the bias point. For some characteristic parameter values 
shown in Table [l] we find that equation ^ gives a phase-separation length i = 120 /im. 

We can also use the simulation to obtain an estimate of the phase-separation length. To do this, we 
run a simulation of 100 TESs in parallel at a constant bias current, and count the number of channels with 
phase-separated TESs. Finally, since the phase separation depends on both and the length of the TES, 
we can use the simulation to map out regions where the TESs are phase-separated or phase-uniform. This 
is shown in Figure |9] The result is in general agreement with equation ([3|. It is important to note, however, 
that Figure |9] assumes a fixed value of ATc for all values of Tc and TES lengths, so each point on the curve 
corresponds to a different value of a. 

6. Effect on Phonon Simulation 

The TES simulation is part of a broader Monte Carlo simulation of the SuperCDMS detectors, which 
models the full charge and phonon physics of each event [HI [18]. Phonons are parsed into individual TES, 
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Figure 9: Simulated phase-separation boundary. The shaded points above the Une are combinations of TES length and Tc that 
produce phase-separation when biased at 1/3 of the normal resistance. Points below the curve produce phase-uniform TES. 



from which pulses are constructed. Modeling the TES using correct parameters is important for producing 
realistic phonon pulses. The effect of this tuning procedure is depicted in Figures [10(a)] and [lO(b)| which 
show some phonon pulses from data and the best simulated matches from simulated template pulses. Before 
the tuning of TES parameters, there was significant disagreement between data and simulation at short 
timescales. After the correct parameters are used, the agreement is considerably improved. 



7. Conclusion 

We have described a simple numerical simulation of phase separation in TESs as used in the SuperCDMS 
experiment. The simulation reproduces the general behavior of Ib-Is curves for ^ 500-element TES arrays. 
The simulation also correctly reproduces the current response to square wave bias voltage, and it can be used 
to estimate the phase separation length at various critical temperatures and TES lengths. This agreement 
is useful for tuning the simulation parameters. 

This simulation has been integrated into the broader Monte Carlo simulation of SuperCDMS detectors. 
Proper TES modeling has been shown to be important for correctly reproducing the behavior of phonon 
pulses on short timescales. More broadly, the TES simulation is an important tool for optimizing the design 
future TESs and detectors. 
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(a) Comparison of simulated (dashed red) and data (solid blue) phonon pulses in a 
SuperCDMS detector, using the tuned TES parameters. Different panes correspond 
to different channels of the detector. 
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(b) The first 400 ^s of each phonon pulse in (a), showing good agreement in the rising 
edge. 



Figure 10: Phonon pulses in data and simulation. 
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